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Abstract 

The vev's of the magnetic order-disorder operators in QCD are found with 
an explicit calculation using the first order formulation of Yang-Mills theory. 



INTRODUCTION 

The definition of non local observables plays a fundamental role in the study of 
YM theory and QCD vacuum. Indeed the well known Wilson line operator gives one of 
the most widely used confinement criteria, namely the area law behaviour of its vev, 
associated with a linearly rising confining potential between static probe charges in the 
QCD vacuum. 

The structure of the confining vacuum has been described by means of the con- 
densation of "magnetically" charged degrees of freedom [|I|, ^ |3], giving rise to the 
so-called dual superconductor model. Several hypothesis on the nature of these config- 
urations and on the dynamical mechanism which leads to their condensation have been 
formulated, and even if this picture in not yet uniquely determined they are nevertheless 
believed to play a major role in the phenomenon of confinement. 

Actually this structure of the QCD vacuum admits different possible phases ||, 
not all confining, which can be labeled by means of the vev's of the Wilson loop W and 
of an other non local operator M, the 't Hooft magnetic disorder parameter. W and M 
give rise to the well known 't Hooft algebra |5|], derived by the implicit definition of M 
as producing center valued singular gauge transformations along the magnetic vortex 
lines. 

The picture of the superconductor model for QCD is analytically realized in the 
first order formulation of YM theory, where an explicit definition of the color magnetic 
operator M is easily given and the calculation of the vev of both M and W can be 
performed, displaying the expected behaviour for the confining phase. 
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The first order form of pure Yang-Mills is described by the action functional [0, [|, [j| 
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Sbf-ym = J Tr(?5 AF + |BA*B) (1.1) 

where F = \F^ v T a dx^ A dx v = dA + i[A, A], D = d + i[A, ■} and B is a Lie valued 
2-form. The generators of the SU(N) Lie algebra in the fundamental representation 
are normalized as TrT a T 6 = 5 a b/2 and the * product (Hodge duality) for a p form 
in d dimension is defined as * = e ll '" td /(d — p)\. The classical gauge invariance of 
( |1.1| ) is given by 5 A = DA Q , SB = — i[Aq, B] and the standard YM action is recovered 
performing path integration over B or using the field equations 

F= l -i-*B , DB = . (1.2) 

Off shell B does not satisfy a Bianchi identity and this fact should be related with the 
introduction in the theory of magnetic vortex lines. We remark that the short distance 
quantum behaviour of ( |1.1| ) is the same as in standard YM as it has been explicitely 



checked [12 



The action functional defined by the first term in the r.h.s. of (|1 . 1| ) is known as 
the 4D pure bosonic BF-theory and defines a topological quantum field theory JTT] . 
Then the bosonic YM theory can be viewed as a perturbative expansion in the coupling 
g around the topological pure BF theory, the explicit topological symmetry breaking 
term ~ g 2 B 2 in (|1 . 1| ) introducing local degrees of freedom in the topological theory. 

The presence of the Lie-algebra valued two-form B field in ( |1 . 1| ) allows the defini- 
tion of an observable gauge invariant operator 

M(E,C) = Tcexpiik Jj 2 y Hol|( 7 )£(y)Hol*( 7 ')} , (1.3) 

where H0IK7) denotes the usual holonomy along the open path 7 = 7^ with initial 
and final point x and y respectively, Hol|( 7 ) = Pexp(i J? dx^A^x)), where k is an 
arbitrary expansion parameter, x is a fixed point over the orientable surface E G M 
and the relation between the assigned paths 7, 7' over E and the closed contour C is 
the following: C starts from the fixed point x, connects a point y G X by the open path 
7s y and then returns back to the neighborhood of x by j' yS , (which is not restricted to 
coincide with the inverse (7m,) -1 = Jyx)- From the neighborhood of x the path starts 
again to connect another point y' G E. Then it returns back to the neighborhood of x 
and so on until all points on E are connected. The path C s = {7U7'} is generic and we 
do not require any particular ordering prescription as it is done in similar constructions 
devoted to obtain a non abelian Stokes theorem |T3] . Of course the quantity ( |1.3| ) is 
path dependent and our strategy is to regard it as a loop variable once the surface E 
is given. 

Using the hamiltonian formalism it is possible to show that M(C) generates 
a local singular (or equivalently a multivalued regular) gauge transformation, fic(x), 
along C: this is precisely the defining property of the 't Hooft color magnetic variable 
[f|. Using the classical constraints which arise from the action (|1 . 1|) it is possible to 
generate the classical gauge transformations; in order to have first class constraints the 
field content of the theory has to be enlarged, including a Lie valued vector field 77. 
This field corresponds to "topological" degrees of freedom which in our case become 
dynamical ||14||. 



Given the classical algebra of tranformations, when switching to operator valued 
quantities, and considering the ordering procedures required by quantization,, one ob- 
tains (assuming S ~ S 2 and k small) 

M(C)\A >~ Tr{E + 2ik <f dy l f da r ^e irs 5 {3 \y - x)5 ab T a T b }\A > . (1.4) 



The integral in (|1.4| ) can be expressed as 



/ dy* [ da r ^e irs 5^(y-x) = -i dy l 



s 4.71 trs dx l 







dxi \ \y — x\ 



~ lim-^ / dy* I dx r e ir i^-( -J—) = Link(C, C') e ^ = sLink(C) , (1.5) 
e-toAnJc Jc oxi\\y — x\J 

where e gives a point splitting regularization between C and C, where S = £' U 
X", <9X' = C and C" encloses the singularity of (|1.5|) . The above linking number is 
the so-called self-linking number of C |TS| . While the linking number of two separated 
loops is not an invariant quantity in 4D, the quantity in (|1.5| ) is well defined and finite, 
takes integer values and equals the number of windings of C around C. 
Putting the above formulae in ( |1.4j) , we find 

M(C)\A >~ Tr{E(l + 2ikc 2 (t)sLink(C))}\A > , (1.6) 

where S ab T a Tf, = C2(t)l. Eq. ( |1.6|) implies that M{C) generates an infinitesimal multi- 
valued gauge transformation; whenever C winds n = sLink(C) times around C, M(C) 
creates a magnetic flux || 
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The finite multivalued gauge transformation fict^] generated by the action of M(C) 
over some state functional is given by 

M{C)\A{x) >= ^^[xK^f) +idg)fl c [x} >~ Tr{e i9 * c E}|A(£) > . (1.8) 

Because of the multi- valued nature of since A Uc = Q^i^A + d)VL c should always 

be single valued, Qq must be in the center of SU(N). To recover the standard form of 
the center, we normalize the free expansion parameter as k = 27r/Nc2(t). With these 
normalizations the form of the color magnetic flux is given by $c = 27m/ Ng and the 
't Hooft algebra is easily recovered || . 



COMPUTATION OF < M(E, C) > 

In this section we compute the average of the BF-observable M(E, C) and precisely 
we consider the normalized connected expectation value < M(C) > conn = <A 5^ > . 

In order to perform calculations we assume the scheme of the abelian projection 
gauge, in which SU(N) is partially gauge fixed to an abelian subgroup 0. In general 
this should be implemented using an interpolating gauge |16| . Choosing (for example) 
Yij = (Fi 2 )ij/(\i — Xj), i ^ j — 1, . . . , N, where the A's are the eigenvalues of F, the 
proposed gauge condition is Y ch + £D° * A ch = where the superscripts ch, stand for 



the off-diagonal and diagonal part of the matrix A = A + A ch and D° is the covariant 
derivative with respect to the diagonal part of the gauge field. Interpolating gauges 
are such that for short distances the relevant gauge is D° * A ch = and the theory 
is renormalizable, while for large distances the gauge is Y ch = 0. In standard YM 
theory we can find in the adjoint representation only a composite of the F^ v (and its 
covariant derivatives), but the dependence of and D° * A ch from the momentum 
is the same. Therefore the dominance argument doesn't apply. This is why the field Y^ 
was introduced. Unfortunately this fact makes this gauge difficult to implement using 
standard quantum YM fields. 

Quite remarkably, these problems are not present in the BF-YM theory due to 
the presence of the microscopic B field. The interpolating gauge can now be easily 
implemented choosing, for instance, B^ + -D * A ch = 0. Equivalently in our formalism 
we can diagonalize the two-form B on the surface S, B = V~ X BV = diag(f3i, ...,/3n), 
/3i > $2 > ■ ■ ■ > /3jV; with V being the singular gauge transformation needed to imple- 
ment the abelian gauge, and then use the background gauge condition D° * A ch = 
in the renormalization program. Indeed this is the way in which magnetic config- 
urations, related to the abelian unbroken group, should enter the theory in the large 
distance regime and fill the vacuum, while in the same limit we may neglect the massive 
off-diagonal degrees of freedom B ch , A ch . This approximation is often called Abelian 
dominance || [TO . 



The existence of monopoles in the Abelian projection gauge is due to the compact- 
ness of the U(1) N ~ 1 group and it is related to the existence of non trivial topological 
objects for the entire SU(N) theory. The reducibility of the SU{2) gauge connection, 
implies that the gauge bundle is split and thus requires the existence of a positive defi- 
nite first Pontrjagin class and intersection number (which we will define later on). This 
fact, in its turn, implies the absence of anti-self-dual harmonic (closed) two forms |T7 . 



In the general SU(N) case we now rewrite the functional integral in terms of the 
variables = [v4 ]jj, i — 1, ...,N, J2i a i — and = [F°]a = d[A°]a. It is convenient 
to rescale the previous geometrical fields A , B° to the physical ones A — > gA°, 
B° — > ^B°. Furthermore, we replace the surface integral in the definition of M with 
the integration over the so-called Poincare dual form, uo^-, of (the homology class of ) 
the surface S 0. By definition, u^, is closed. Moreover, choosing an orientation of 
the four manifold and remembering the absence of anti-self-dual harmonic two forms 
(imposed by requiring the existence of a reducible gauge connection) u; s is chosen to 
be self-dual i.e. = uJs and with the property that (up to gauge one-forms) for a 
generic two-form t 

f t~ I uj^At . (2.1) 

In a local system of coordinates (x,y,u,v) on E, so that S is given by the equations 
x — y — 0, the dual form u; s can be taken as ujy, — 5^(x, y)dx A dy and normalized as 



/ w s 

Jn{e) 



5 {2 \x,y)dxAdy = l , (2.2) 



where iV(£) is the transversal tubular neighbourhood on the surface S. 

Using the Abelian dominance approximation and defining the expansion parameter 
k in units of the bare color charge g with a suitable normalization, k = 2nqg, the 't 
Hooft loop operator M(C) becomes 

M(C) = Tr{O lJ (C)} , (2.3) 



A 



Oij(C) = Sij exp{i2nq / w s A fa[cos(g f aj) + i sin(g f aj)]} , 

JM 4 JC JC 

where the configurations ctj are to be intended as singular ones, related to the singu- 
larities occurring if two consecutive eigenvalues of B coincide [0]. If fa = fa+i, we shall 
label such a point by and qi will be the corresponding magnetic charge. The key 
point will be the identification of the arbitrary expansion parameter q with the magnetic 
charges qi, i.e. we will set qi oc q oc K Consider then the magnetic order parameter 
M in the q — > limit, corresponding therefore at the strong coupling limit g — > oo. In 
this limit, and relaxing the constraint J2i a i — thus extending the maximal torus of 



SU(N) to U(l) Hl8| ) the operator ( |2.3|) can be approximated by 



M(C)=20«(C7)~^exp{^y'u^A2Acos(^jf Oi)} . (2.4) 

With all these approximations taken into account the vev of the magnetic order 
parameter in the strong coupling region becomes 

N r 1 . /" 1 47TO 

< M(C) > conn ~ — ^ y exp{-- J2 J A A *ft + zft A [doj + -^w s ■ 

■cos(#^ a*)])} = < ^ > y £>a;exp{--^ J [dojj + cos(# ^£ c^)] 2 )} , (2.5) 

where the square of a form t means t A *t. We now split the gauge fields as ai = 
cti + hQi, where the quantum fluctuations Qi must be gauged (e.g. by a covariant gauge 
condition) and the a^ are singular classical configurations. Postponing for a while the 
discussion of quantum fluctuations, we concentrate on the semi-classical contribution 
to the path integral which is given by 



< M(C) > conri ~ iVexp{ — j^-J2j u ^ Awjycos^j^OiJcos^^aJ)} , (2.6) 

where E' represents a point splitting regularizations of E for any point of E. The above 
equation has been obtained observing that the configurations a« obey the (monopole) 
equation 

* don = -irr^T, cos(# * , (2.7) 
JM Jc 

derived using the self-duality and closedness of u^. Due to the singular behaviour of 
the ctj's partial integration is not allowed |TJ| and no electric currents are present in 
the model. Equations of the type (|2.7|) appeared already in the study of the duality 



properties of gauge theories and 4D manifold invariants |21 
The classical contribution to the action is 



27r 2 ?7i 2 f , f r 2,7r^Tn? N 

S = — / J! A u-z> cos(# j> ai) cos(# j> a[) = Q(E, E') , (2.1 

g J a Jc jc g 



where Dirac quantization g § c 6ti = 2im has been used, m G Z labels monopole charges 
and Q(H, E') is the algebraic intersection number The closed curve C = {y M (£)} 
corresponds to the framing contour of C = {x M (t)}, i.e. if it happens that y^(t) = 
x^lt) + en M (t), with e — > 0, and = 1 where n^(t) is a vector field orthogonal to 



c 



C. Then <3(E, E') becomes the self-linking number of C, Q(E, E') = sLink(C). From 
a physical point of view, we may define sLink(C) as 



sLink(C) 



L(C) 
P 



(2.9) 



where L(C) is the perimeter of the loop C and p plays the role of unit of lenght, 
reasonably associated to the vortex penetration depth. 

Let us now discuss quantum fluctuations. If the effective theory for large distances 
is a U(l) type theory, for short distances the charged degrees of freedom cannot be 
discarded anymore. Let A be the scale separating these two regimes and let us divide 
the gauge field accordingly; moreover let, for the sake of simplicity, the gauge group 
be SU(2). For scales bigger than A we take A 3 = a + Q 3 which is the usual U(l) 
prescription. For scales smaller than A we take A a = 5 a3 a + Q a , i.e. we continue the 
classical solution into the small scales region where the quantum fluctuations coming 
from the charged degrees of freedom cannot be discarded. The expectation is that the 
small scales behaviour is insensitive to the classical solution according to the background 
field method. Performing the functional integration over the quantum fluctuations leads 
to a double contribution, in complete analogy with the saddle point evaluation around 
an instanton background ||23|| . The first contribution is given by a ratio of determinants 
given by 



R 



Def(-L 



where D = d 



iga. 



D 2 



Det'(- 



-L) 



Det(-D 2 ) 
Det(-d 2 ) 



;2.10) 



D 2 5„ 



[1 - l j)D,D u + ligFl 



and L is 



given from L evaluated around the trivial background. £ is the gauge parameter and 
in the primed determinants zero modes are omitted. 

The second contribution is given by the Pauli-Villars regularization of the determi- 
nants and it amounts to a scale \x (which is the Pauli-Villars mass) raised to a certain 
power which is given by the dimension of the moduli space of the classical solution. 

Let us now proceed with the evaluation of these two contributions. Using the self- 
duality property of our classical solution ( |2.7|) , the ratio of determinants Q2.10 ) can be 
written as R = \Det(-d 2 ) / Det(-D 2 )] @. This ratio has been evaluated in Ref.[g§] 
using the heat kernel method in the case of an SU(N) gauge group but it is easy to 
specialize this result to our case. We obtain 



In 



Det(-d 2 ) _ ln(/iA) 
Det(-D 2 ) ~ 96tt 2 



~*L(x)fI,(x)#x = X 4^r- A I (^u*)(cfcU*K* . 



96vr 2 4 



(2.11) 

The factor 1/4 comes from the normalization of the gauge group generators according 
to the reduced connections «j. 

The contribution coming from the regularization of the zero modes is obtained 



once the dimension of the moduli space is computed, according to Ref.[|2T[], to be 



dim M = ci(L) 2 = ci(L) A c x {L) 



32tt 2 



(da) ^ (x) (da) „ v (x) d 4 x 



(2.12) 



Putting together the classical result ( |2.8| ) with the quantum fluctuations, we find 
that the bare coupling g can be substituted by its renormalized expression and that 



the exponent in ( |2.6| ) can be written as 



8vr 2 Cl (L) 2 Cl (L) 2 /8vr 2 2 \ _ Cl (L)~ / ^tt 



> , , , , > =^-(^-/?iln(/xA)) , (2-13) 

0a 8 8 V # 2 3 / 8 V # 2 / 

where ft = y is the first coefficient of the SU(2) beta function of the non-abelian 
Yang-Mills theory. 



THE AVERAGE OF THE WILSON LOOP 

In this section we shall compute the average of the Wilson loop and find an area 
law behaviour for its leading part. Furthermore, in our formalism, the area law gets a 
nice geometrical interpretation: it is the response of the true QCD vacuum to arbitrary 
deformations of the quark loop C. 

The starting point here is given by the Wilson loop operator written in terms of 
the non abelian Stokes theorem (see e.g. 



W t (C) = W t (Z,C) = Tr t {P s exp[z jf Hol|( 7 )F(x)Hol*( 7 ')]} , (3.1) 

where C = <9£, C = U j'( x )} was defined at the beginning of section 2 and 

Ps means surface path ordering. W t is calculated with respect to some irreducible 
representation t of SU(N). 

<W t (C)> conn = — = JVBVA ~ Sbf—ym ' ^ 2 > 

where Sbf-ym was defined in Expanding in series the Wilson loop (|3.1j) we get 

< W t (E,C) >= E < V r * Ps / • • • / (^Hol( 7 )F(a;)Hor 1 (7 / )) n > ■ (3.3) 

We then use the identity 

e -U* B ^FUx)=4i S ——( e -U , (3.4) 

6*B a pa (x) 

and performing a partial integration with respect to the functional derivative in (|3.4j) 
we can replace, in the path integral, 

^Hol( 7 )e-i/* i?Q - Fa F(x)Hor 1 (7 / ) -> 4Hol( 7 )e-^ / * g "- Fn ( = - , J HoP^V) . (3.5) 

* -B(x) 

The functional derivative acts only to its right on the exponential of the mass term 
— f? 2 / 16 / B^ U B^ U , since Hoi does not contain the B field. We now need the identity 

where V is the functional defined by 



Resumming the exponential series for the Wilson loop we finally get the "duality" 
relation 

<W l( C)>^ < M ^f 1 C -^> , (3.8) 

where 

2 

M t *(E,C,C = <9E) = Tr,[P s exp{-|- jf Hol5( 7 ) * B(aOHo£(Y)}] • (3.9) 

M t *(E,C) is the dual (in the sense that B — > *£?) of the observable M t (T,,C) defined 
in flOD with k set to k = ig 2 /2. 

To calculate Q3.8| ) we expand perturbatively in (7 both the exponential and the 
holonomies which appear in the exponent of (|3.9| ) p5|| . The first relevant contraction 
encountered at lower level is given in terms of < A * B >, which can be computed 
starting from the off-diagonal propagator < AB > fl2 |. Therefore we find 

<M* t {C) > conn =e-^ C2W ^/ s <^>A(S) . (3.10) 

A(E) depends on higher order integrations over E which do not simply involve the 
quantity § c ^ < A * B > alone; therefore while the explicit calculation of A(E) is an 
open problem, for the purpose of showing the area law behaviour its knowledge should 
not be essential. Consider now the integral in the exponent of ( |3.10| ), 

[ <A*B >=< I A f da^(x)(*B(x))au > (3.11) 
c jt Jc Jt 

A 



=<f A (*dar(x)B(x) flv >= i dy A / (*da)^(x) < A(y) x B(x))^ > , 

where *da(x) is the infinitesimal surface element of the plane E* dual to E at the point 
x G E. We may rewrite 

I dy\*daY\x) < A(y) x B(x)) flv >= I dy x [ < A X B^ > , (3.12) 

recalling that uj^, is locally given by a bump function with support on E. Eq. ( |3.12| ) is 
by definition the linking number between the curve C and the dual plan E* in x to E. 
Indeed the linking Link(C, E), whith arbitrary C and E, is defined by |26 



Link(C, E) = -L / dx a [ da^(y)e 
8tt^ Jc Jt 



a/375 1 1 

TT* JC JT, \x — y\ 



5 



= A^dx a Jja^(y)<A a a (x)B a py (y)> . (3.13) 

In our case, by construction, Link(C, E*) 7^ 0. The residual integration over E in ( |3.12| ) 
spans all the dual E* to E, yielding a contribution proportional to the area of E, 

<W t (C)> conn ~ exp{-^^ /Link(C,E:)} , 

o JT 

Link(C,E*.) ~ A(E) . (3.14) 



We may get a better understanding of ( p,14| ) considering a lattice regularization of E, 
i.e. E — > Ep£. In this case C runs over the links of the lattice, while E* corresponds 



to an element of the dual lattice through the plaquette centered at x. J EpL Link(C, £*) 
is an integer which counts the number N v of vertices of the dual lattice on Ep^ or 
equivalently the number of plaquettes Np of Ep^. We may then write 

/ Link(C,E*) =N P = ^p^- , (3.15) 

where A(Ep^) is the minimal area bounded by the quark loop C and a is the lattice 
spacing. Up to this point the quantities which enter the calculations are the bare 
ones, but when passing to the continuum limit a — > the renormalization of the theory 
implies that they are replaced by the renormalized ones. In particular Np = AiT^pp) / a 2 
should become A(E)/l 2 with I a typical scale in QCD which, owing to a proper dressing 
of the off diagonal propagator entering fl3.10|) , should be choosen as Aq CD . Therefore 
one may rewrite ( p,14j ) as 

< W t (C) > conn ~ exp[-a(A QC p)A(S)] , (3.16) 

where <t(Aqcd) is the string tension defined by 

a(A QCD ) = 9 -^A 2 QCD (3.17) 

and where we have replaced the bare coupling constant g with the running coupling 
constant gp. The way in which the continuum limit is reached and the fields are dressed 
is presently not under our control but the key geometrical feature of the area law should 
remain unchanged. 

A deeper investigation should also clarify some other questions. While the area law 
displayed in fl3.16| ) corresponds to the expected confining behaviour of QCD vacuum, 
no screening is shown when the test charges are taken in the adjoint representation and 



the behaviour in this case is closer to what expected in the large N limit [27fl. It might 
seem puzzling that the area law behaviour comes out of a perturbative calculation. 
In reality our perturbative expansion has to be understood as an expansion around 
the gauge fixed topological BF theory whose vacuum contains "relevant" topologically 
non-trivial configurations. It is thus resonable to imagine that the non perturbative 
information comes out of this non-trivial vacuum structure. The contribution of these 
"magnetically charged" configurations is expected to give rise to a mass gap in the 
model and we guess a possible form for it, namely 

m 2 ~ fi 2 e~^ , (3.18) 

where \i is a regulator mass and q is the magnetic charge satisfying a proper Dirac 
quantization qg = 27m. A suitable dressing displaying such a pole for propagators 
is obtained substituting (in euclidean momentum space) the 1/p 2 denominators with 
l/(p 2 + U(p 2 )), where IT(p 2 ) ~ 2(|^)p 2 ln(^). In such a way, in a q —>■ limit (i.e. 
strong coupling expansion) and applying Dirac quantization, the string tension read 



out of (|3.14T) should be vanishing when the adjoint representation is choosen. 

A last point to be discussed is the apparent surface dependence which the Wilson 
loop acquires in the representation (|3.8|) . Deforming the surface bordered by the loop C 
corresponds to a "gauge" transformation on the connections defined on the loop space; 
in this sense the "gauge" independence of observables in loop space corresponds to 



n 



surface independence for the Wilson line. Lastly, our formalism appears quite similar 
to that recently introduced in [^8| for the vev of W. 
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